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Abstract 

We report a new model for highly mismatched semiconductor (HMS) alloys. Based on the 
Anderson impurity Hamiltonian, the model generalizes the recent band anti-crossing (BAC) model, 
which successfully explains the band bowing in highly mismatched semiconductors. Our model is 
formulated in empirical tight-binding (ETB) theory and uses the so called sp 3 s* parameterization. 
It does not need extra parameters other than bulk ones. The model has been applied to BeSe^Tei-^ 
alloy. BeTe and BeSe are wide-band gap and highly mismatched semiconductors. Calculations show 
large band bowing, larger on the Se rich side than on the Te rich side. Linear interpolation is used 
for an arbitrary concentration x. The results are applied to calculation of electronic and optical 
properties of BeSeo.41Teo.59 lattice matched to Si in a superlattice configuration. 

PACS numbers: 73.21.Cd, 71.55.Gs, 78.66.Hf 
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I. INTRODUCTION 



In recent years there has been a tremendous interest in electro-optical properties of II- 
VI semiconductors in particular epitaxial II- VI heterostructures. New developments^ in 
the growth of Si lattice- matched BeSe .4iTe .59 open the opportunity to a new class of Si 
based devices. Be-chalcogenides are wide-band gap zinc blende semiconductors with lattice 
constants close to that of Si. Thus BeTe and BeSe have the lattice constants of 5.6269 
and 5.1477 A, respectively, 3.6 % larger and 5.2 % smaller than Si. Vegard's law says that 
the lattice matched composition with Si is BeSe .4iTe .59. Therefore Be-chalcogenides are 
candidates for Si-based heterostructures. 

The difference in size and orbital energies between Se and Te in addition to large lattice 
mismatch between BeTe and BeSe makes the virtual- crystal approximation (VCA) inappro- 
priate for the ternary alloy BeSe x Tei_a;. The band anti-crossing (BAC) model^ has been 
introduced in order to explain the electronic structure of highly mismatched III-V-N alloys 
and II- VI alloys like ZnSe x Tei_a^. At impurity like concentrations close to both end points, 
the electronegativity difference between constituent elements gives rise to localized energy 
levels close to the conduction or valence band. Thus in the ZnSe(ZnTe)-rich side, the band 
gap bowing is mostly determined by the anticrossing interaction between the Te(Se) local- 
ized level, which behaves like an impurity, and the extended states of ZnSe valence band 
(ZnTe conduction band) near the center of the Brillouin zone^. 

In this communication we develop a model which is a natural extension of the BAC 
model to empirical tight-binding (ETB). Based on this model we determine optical bandgap 
of BeSe x Tei_ x alloy and further we analyze the band folding in Si/BeSe .4iTe .59 heterostruc- 
tures. 

II. MODEL 

There are several studies employing tight-binding (TB) models for HMS. They use either 
supercells 6 or add extra parameters to the usual TB parameters within the BAC modet- Our 
approach needs no extra parameters others than the usual TB parameters and is a natural 
extension of the BAC model to TB. We consider first the dilute limit. The starting point is 
the impurity model of Anderson^. Consider a complex-structured impurity interacting with 
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the host crystal having the Hamiltonian 



H = H C + H imp + V = H + V (1) 

where H c is Hamiltonian of the host crystal, H imp is the Hamiltonian of the impurity and V 
is the Hamiltonian of the interaction between impurity and crystal. The above Hamiltonians 
have the following expressions 

Himp — ^ ] £da C pda C pdcn (2) 
p,d,cr 

V= {^kdC+dvdka + h.c). (4) 

p,d,i,k,a 

Here a is the spin index, d are the energy levels of impurity at the pth site (N iy the total 
number of impurities), i is band index of the crystal, k is the wavevector, Cp da {c p do)axe the 
creation (annihilation) operators of electrons on impurity levels, c^^c^) are the creation 
(annihilation) operators of electrons of the bands in the crystal, and V is the coupling 
between impurity and crystal and has factor 1 j \f~N with N total number of sites. In order 
to calculate the effect of impurity we shall use an expression for the projection of H into H, 
a subspace of the Hilbert space spanned by the Hamiltonian H . We define P (and Q) the 
projections onto (out of) H as P 2 = P, Q 2 = Q, P + Q = 1, PQ = 0, [P, H ) = [Q, H ) = 0. 
Let the Green function of the entire system be denoted by G (z) = (z — H)~ and the 
projected part onto H as G (z) = PG (z) P. In this way we denote the projection onto H of 
any operator A as: A = PAP. Straightforward algebra gives us 

G(z)= 1 , (5) 

z — H o — R (z) 

where R (z) = PVP + PV z _ h ®_q V q VP. One can expand R (z) for small V 
— , s Q Q Q 

R(z)=PVP + PV ^ TT VP + PV ^—V VP + 6 

z — H z — H o z — H Q 

Due to Q in the expansion, all the intermediate states are outside of H, therefore the 
diagrams representing R (z) must be irreducible. If R (z) is small compared to H c , Eq. (jSJ) 
can be expanded in a power series of R (z) as 
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G (z) = G C + G C RG C + G C RG C RG C + .... (7) 

No approximations have been made so far. If V is small we can replace R (z) in Eq. 
([6]) by the first two terms and summing up all contributions (the first one is in our 
model). By making such an approximation we, in fact, sum an entire class of diagrams. 
The impurity averaging is made by noticing that all macroscopically observables are self- 
averaging, i.e. they have asymptotically exact values in the thermodynamic limits. In other 
words the average of the product of such quantities is equal (within asymptotic accuracy) 
to the product of their averages. Therefore the impurity averaging in Eq. (jSJ) is simply 



taken as averaging R(z) in Eq. flBj). By keeping the first two terms in Eq. (j5J), the model 



is equivalent to the optical model laid out in Ref. LLO: it only brings the shift to the energy 
levels of the unperturbed Hamiltonian and gives the following form for the diagonal part of 
the Green function of ith band 

Gkkaa ( Z ) = ~Z , v ,2 ) (8) 

1 _ _ _ I Km I 

1 tika X z _ £d<T 
a 

where x is the dilute concentration. For one band and 1-level impurity the result is iden- 
tical to the BAC model 3 -^^. This result can be easily expanded to include pair impurity 
interactions in addition to single impurity interactions^. 



III. APPLICATION TO SP 3 S* TB HAMILTONIAN AND NUMERICAL RE- 
SULTS 

The model laid out in the preceding section is directly applied to sp 3 s* Hamiltonian^ 3 - 
with spin-orbit interaction^ The TB Hamiltonian is written in the sp 3 hybrid basis and the 
basis is rotated in such a way that a unit cell is formed by the anion hybrid orbitals and the 
cation hybrid orbitals pointed toward the anion site. The s* orbitals remain unchanged. In 
the no spin-orbit case, the transformation is Hh y b (k) = S + (k) Hao {k) S (k), where Hao (k) 
is the TB Hamiltonian in atomic (Lodwin) orbital basis and (k) the Hamiltonian in 
sp 3 hybrid basis. The S matrix is block diagonal in anion/cation index and has the form for 
anion and cation site, respectively, as 
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2 



11110 
1 1 -1-10 
1-1 1 -10 
1-1-1 1 
1 



(9) 



1 
2 



-ik- d\ 



-ik- d,2 



^,—ik-di Q—ik-d2 



_g— ik-di g— ik- d2 



-ik- d'j 



-ik- d$ 



e ~ik-d 4 g 
e ~ik-d 4 g 



_g— ik-dz g— ik- d 4 g 



_ e ~ik-d 4 g 

1 



(10) 





The vectors in (TTOT) are (a is the lattice constant): di = (a/4) (1, 1, l),d 2 = (a/4) (l, 1, l) ,d 3 = 
(a/4) (l, 1, l), and d 4 = (a/4) (l, 1, l). An impurity is such a cell interacting through hybrid 
orbitals with the average/host crystal. 

BeSe and BeTe are quite new materials. The nearest-neighbor sp 3 s* parameters were 
fitted to GW calculations 15 . They reproduce valence band edges r 8 and T 7 , and the conduc- 
tion band edges Tq and X\. The average crystal is considered by the average parameters. 
The hopping parameters were scaled according to the Harrison scaling rule^ and then av- 
eraged. The band offset between BeTe and BeSe is considered to be 0.41 eV as indicated 
in RefJl7l. In this way one calculates the bowing of the on-site energies in addition to linear 



terms given by VCA. The hybrid states of the impurity lay outside bandgap of the host 
crystal, such that the net effect is large deviation from linearity of the band edges of the 
alloy. Mathematically one diagonalizes an extended Hamiltonian and accounts for the shifts 
in the band edges at the V and X point in the Brillouin zone. Those shifts are used to 
calculate the bowing parameters for each self-energy. 

We calculate the bowing parameters of the direct and indirect bandgap around and 
1 limits of concentration according to BAC model and follow the spirit of the VCA to 
interpolate linearly the effect of BAC model between these limits. The linear interpolation 
has been successfully used to fit experimental data for ZnSeTe alloy^. Linear interpolation 
for the direct and indirect bandgaps is consistent with linear interpolation performed on 
the tight-binding parameters. The tight-binding parameters for BeTe and BeSe are shown 
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TABLE I: Matrix elements in eV of nearest neighbor sp 3 s* model including spin-orbit interaction 



for BeTe and BeSe. The on-site energies of BeTe have b een upgrad ed by 0.41 eV, the band offset 



Vogl et al^ . The bowing parameters of 



between BeTe and BeSe. The notation is according to 
the on-site energies for BeTe-rich limit and BeSe-rich limit are shown in parenthesis. For a general 
concentration x the linear interpolation is used. 





BeTe 


BeSe 


E(s c) 


5.112+0.41 (-1.85) 


5.560 (0.55) 


i-j y o . ex i 


_i 5 4m _i_n 41 (0 R4) 


-14 953 f1 




4 427+0 41 (-Q 6) 


5 026 f-5 7) 


EfD a) 


_n 299+0 41 (Q 5) 


300 (5 81 


E(s* c) 


30.16+0.41 (1.0) 


21.666 (1.3) 


Efs* a) 


39 203+0 41 fO 5) 


24 433 fO 65) 


V(s,s) 


-3.303 


-8.195 


V(sc,pa) 


4.423 


5.633 


V(sa,pc) 


5.511 


4.89 


V(x,x) 


0.331 


1.531 


V(x,y) 


6.362 


6.324 


V(s*a,pc) 


11.503 


7.462 


V(s*c,pa) 


3.11 


4.572 


A a 


0.97 (-0.4) 


0.499 (-0.15) 


A c 









in Table [H The bowing parameters of the on-site energies for Te-rich and Be-rich limit, 
respectively, are shown in parenthesis. 

Calculated direct and indirect bandgaps are shown in Fig. Q] for VCA and BAC models 
against linear interpolation. The conduction band minimum is located at the X-point in the 
Brillouin zone, such that the fundamental bandgap is indirect. The VCA model gives almost 
constant bowing parameters of 0.49 eV for the direct (optical) bandgap. The bowing of the 
direct bandgap given by the BAC model is much steeper on the Se-rich (9.8 eV) side than on 
the Te-rich side (2 eV) suggesting that utilizing just one bowing parameter is inappropriate 
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FIG. 1: Direct and indirect bandgaps in BeSe x Tei_ x : linear interpolation (full line); VCA (line 
with squares); B AC (diamonds); and interpolation to BAC model (dotted line). 

to describe bandgaps in these structures. The results for the indirect bandgap follow the 
same trend as that of the direct bandgap (Fig. [T]) with the minimum of the indirect bandgap 
at 1.7 eV for x around 0.6. This trend is similar because of the large bowing of the valence 
band edge in addition to conduction band edge. Moreover, the VCA results are very close 
to the linear bandgap. 

We calculated also electronic and optical properties of a Si/BeSe .4iTe .59 superlattice 
(SL) in (001) direction. Abrupt interface, flat band conditions were assumed. We adjusted 
the conduction band offset between BeSe .4iTe .59 and Si at 1.2 eV as determined from 
electrical measurements^. Two interface subbands were found, one empty and one occupied 
within the Si bandgap. The origin of these interface subbands is due to polar nature of 
the interface or large difference between on-site energies of Si on the one side and Be or 
Se/Te on the other side.— In Fig. [2] we plot the joint density of states for vertical transitions 
in (Si2)4(BeSeo.4iTeo.5g)4 SL. In such structures, due to band folding, the threshold of the 
direct bandgap is lowered for Si from 3.35 eV to 1 eV. The threshold is slightly below the 
fundamental bandgap of Si because of the two interface bands. Moreover, the first peak 
is wider. Si as an indirect bandgap semiconductor has two kinds of confined states in a 
quantum well 19 , one given by the longitudinal valleys with an effective mass of 0.19m and 
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Energy [eV] 

FIG. 2: Joint density of states responsible for vertical transitions in Si/BeSe^Tei-^ SL. 

the other given by transverse valley with an effective mass of 0.91m . Hence the confined 
states that have a different first peak is made of the first states of both types of valley. 

IV. CONCLUSIONS 

In conclusion, we extended the band anticrossing model to the empirical tight-binding 
theory. We used the Anderson model for impurity and a sp 3 hybrid basis for zincblende 
structures within sp 3 s* Hamiltonian. The effective Hamiltonian of the alloy was obtained by 
impurity averaging and keeping only the terms responsible for energy shifts due to alloying. 
The pair impurity effects can easily be included as well. Thus no extra parameters are 
needed to calculate bandgaps. 

The model was used for BeSe^Te^ alloy. BeSe^Te^ shows large band bowing, larger 
on the Se-rich side similar to ZnSeTe system. Bandgap was interpolated linearly between the 
two dilute limits. Further the model was applied to Si/BeSe .4iTe .59 superlattice in (001) 
direction with abrupt interfaces. Due to polarity of the interface, two interface subbands are 
found within bandgap of Si. Also calculations show that the threshold for direct transitions 
is lowered in Si and that the absorption edge is slightly below the Si fundamental bandgap. 
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